Electron Correlation Effects in Hyperfine Interactions in 45 Sc and 89 Y 



1,2 



Bijaya K. Sahoo 1 Thomas Beier, 2 B. P. Das, 2 R. K. Chaudhuri, 3 Debashis Mukherjee 



in 
o 
o 

(N 

c 
3 



1 Atomphysik, 

Gesellschaft fur Schwerionenforschung mbH, Planckstra/3 e 1, 64291 Darmstadt, Germany 
2 Non- Accelerator Particle Physics Group, 
Indian Institute of Astrophysics, Bangalore-34, India 
3 Department of Physical Chemistry, 

Indian Association for Cultivation of Science, Calcutta-700 032, India 
(Dated: Received date; Accepted date) 

The relativistic coupled-cluster theory has been employed to calculate the magnetic dipole and 
electric quadrupole hyperfine structure constants for the stable isotopes 45 Sc and 89 Y. The role of 
electron correlation is found to be very important. The trend exhibited by these effects is rather 
different from previously studied single valence atomic systems. 
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I. INTRODUCTION 

Electron correlation effects play a crucial role in the ac- 
curate determination of hyperfine struture constants . 
We have calculated these constants for several atomic sys- 
tems with single s valence electrons 0, using the rela- 
tivistic coupled-cluster (RCC) theory which takes into ac- 
count the interplay of relativistic and correlation effects. 
Our recent work on Pb + which has a single p valence 
electron £| reveals the peculiar behaviour of the elec- 
tron correlation effects in the hyperfine constant of the 
&p 2 Pz/2 state. It would indeed be instructive to study 
how electron correlation effects influence systems with 
single d valence electrons. In this paper we have carried 
out ah initio calculation of the magnetic dipole (A) and 
electric quadrupole (B) hyperfine structure constants for 
stable Scandium (Sc) and Yttrium (Y) which have 3d 
and Ad valence electrons respectively. Hyperfine struc- 
ture constants for ground and first excited states for these 
two atoms have been measured using different methods 

B nun nana. 



II. THEORY 

The hyperfine interaction Hamiltonian with magnetic 
dipole and electric quadrupole terms can be written as 



H hf = AI.J + B 



3(1.3) 



2 + fl.J-/(/+l)J(J+l) 
21(21 - 1) J (2,1 - 1) 
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where A and B are the hyperfine structure constants and 
I and J are the nuclear and electronic angular momenta 
respectively. In perturbation theory, the first order en- 
ergy shift in the atomic state \JMj) with Mj as az- 
imuthal quantum number corresponding to J due to the 



above interaction Hamiltonian yields 



W F 



AK B 



K(K+1)-2I(I+1)J(J+1) 



(2-2) 



1(21 - 1) J(2J - 1) 
where 

K = F(F + 1) - J(J + 1) - 1(1 + 1) (2.3) 

with 

F = I + J (2.4) 
is the total angular momentum of the system. 

The magnetic dipole hyperfine constant A is given by 



R 



A = hn9i 



<J||T«||J) 



y/j(J +1)(2J+1)' 



(2.5) 



where T (1) = = £ -iey/^ipr^a.Y^ (f) with 

*ka re P resen t s t ne vector spherical harmonic, gi = 
is the nuclear Lande g-factor with [ii is nuclear dipole 
moment and fiN = eh/2m p c. 

The electric quadrupole hyperfine constant (B) is given 

by0 



B = 2eQ 



2J(2J - 1) 



(2J+l)(2J + 2)(2J + 3) 



1/2 



<J||T( 2 >||J(}2,6) 
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where Q is the nuclear quadrupole moment and 
T^=J2t^ =j:-er- 3 C ( q k) (r). 

The reduced matrix elements of the above operators in 
terms of single particle orbitals given by 

(«/ii4 1) ii«*> = -(«/+«i)<-«/ii cr(1) n«<> 

/ dr^(P f Q l + Q f P l ) (2.7) 
Jo r 
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and 



and 



r°° 1 

( Kf \\tW\\ Ki ) = -{K f \\CW\\ Ki ) / dr-iPf^ + QfQt), Ujl^in)) 

Jo r 



1 

n 2 



(2.8) 



a=l 



where i and / represent initial and final orbitals respec- 
tively and 

(Kf\\CW\\ Ki ) = (-l)i/+V2^ (2i/ + 1)(2i . + 1) 

with the condition 7r(Zl, 12, 13) = 1 if 11 + 12 + 13 — even, 
otherwise zero. 



III. METHOD OF CALCULATION : 
RELATIVISTIC COUPLED CLUSTER THEORY 

We have performed our calculations in the Dirac- 
Coulomb (DC) approximation which embodies the 
dominant relativistic and correlation effects. The Breit 
interaction which is about two orders of magnitude 
smaller than the Coulomb interaction has therefore 
been omitted in the present case as the systems are 
comparatively smaller in size. 

The Dirac-Coulomb atomic Hamiltonian is given by 
H = ca. Pj + (J3- l)c 2 + V nuc (rj) + £ — , (3.1) 



3<l 



Til 



where a and [3 are the usual Dirac matrices and V nuc (rj) 
is the potential at the site of the j th electron due to 
the atomic nucleus. The rest mass energy of the elec- 
tron is subtracted from the energy eigen values. We 
first solve the relativistic Hartree-Fock (Dirac-Fock (DF)) 
equations to obtain the single particle orbitals and their 
energies, 

H DF = ca.p 3 + {13- l)c 2 + V nuc ( rj ) + Uj 



. The residual Coulomb interaction is given by 



Til 

3<l 3 3 



(3.2) 



The single particle orbitals are obtained by solving the 
following equation self-consistently 



where 



(tj + Uj)\(f>j) = tjtyj), 



ca ■ pj +(fi- l)c 2 + V nuc (rj) 



(3.3) 



ba(r2)\ — \Mr2))\Mn)) 

Tl2 



occ represents total number of occupied orbitals. 

The single particle relativistic orbitals can be expressed 



as 



u.( r )\ -If P j( r )\XK jmj ) \ 

where Pjir) and Qj(r) are the radial part of the large 
and small components respectively and \xnm ) and 
IX-re-m ) are their respective spin angular momentum 
components, e^'s are the single particle energies. 

We have employed the RCC to incorporate correlation 
effects among electrons due to the residual Coulomb in- 
teraction. In this approach the exact atomic wavefunc- 
tion for the closed-shell system can be expressed as pj 



|*cc) = e T m 



(3.4) 



where T is the core electron excitation operator and 
we call it as closed-shell RCC operator. |$) is the 
above closed-shell determinantal state built out of the 
Dirac-Fock single particle orbitals. 

In the closed-shell coupled-cluster theory one starts 
with the equation 

He T \<f>) = Ee T \^>). (3.5) 

The energy and amplitude determining equations are 

($ K |i7|$) =E8o,k, (3.6) 

where H = e~ T He T , \& K ) is a determinantal state 
with K = 0,1,2.... representing the reference state and 
excited determinantal states. We have considered all 
possible non-linear terms in T- operator for its amplitude 
determining equations. 

Goldstone Q, ItH HE} and angular momentum diagram- 
matic pi 0] techniques are used for evaluating different 
radial integrals and angular factors. The normal ordered 
Hamiltonian is defined as 



H N = H - ($|£r|#> = H - E DF , 
where E DF = 



(3.7) 



We have truncated our wavefunction expansion at the 
level of singles and doubles (CCSD) and all possible 
non-linear terms have been included in the above 
equation. First we evaluate wavefunction for the closed 
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shell system using the above RCC approach and then 
append the corresponding valence electron using the 
open shell RCC (OSCC) method 0,113. 

The new reference state of the open-shell system with 
one valence electron v can be expressed as |l7j 



!<*%)} 



(3.8) 



where a v is the particle creation operator. The ex- 
act atomic states are defined now, using the Fock-space 
OSCC method, as 



I*,} = e T {e s -}|$,}, 



(3.9) 



where S v is the valence excitation operator and we call 
it as open-shell RCC operator. Since systems under con- 
sideration can be treated as with one valence electron, 
the S v - operator exponential series naturally truncates 
at the linear term, i.e. the open-shell wavefunction has 
the form 



=e T {l + S„}|$„>, 



(3.10) 



where 



S v = Sit) + S 2v = ^2 a+a-usP + - a-p a qO-ba v s p V: 



b 



p¥*> 



and 



Sly S dpQ V Sy 



pj^v 



S2v = ^y^ J a^a^a b a v : 



vb 



(3.11) 



bpq 



with and s p v \ are the cluster amplitudes corresponding 
to single and double excitations involving the valence 
electron. 

In the next step, we include approximate triple exci- 
tations by contracting the two-body operator (V es ) and 
the double excitation operators (T2, S^) in the following 
way [3 



qpqr _ 
J vbc ~ 



V es T2 + V es S2v 



e v + e b + e c 



(3.12) 



where q is the orbital energy of the i'th orbital. Note 
that we use notations a,b,c..., p,q,... and i,j,... for 
the core (occupied), particle (unoccupied) and general 
orbitals respectively. 

The equations for the open-shell cluster amplitudes are 
determined from 



($ V \H N {1 + S V }\$ V ) = AE(v) 



(3.13) 



and 



(K\Hn{1 + 5„}|$„) = -AE{v)($* v \{S v }\$ v ), (3.14) 



where AE(v) is the electron attachment energy which is 
equal to the negative of the ionisation potential for the 
valence electron v. 

The expectation value for a general one particle opera- 
tor in a given valence electron (v) state can be expressed 
in coupled-cluster theory as 



(0) v 



(Vv\Q\*v) 
{*v\Q\*v) 

l + N v 

($ v \{l + Sl}e Tt Oe T {l + S v }\* v ) 
l+N v 

($ V \{1 + St}Q{l + S v }\$ v ) 
l + N v 



(3.15) 
(3.16) 



where the normalization term for the v th orbital is ob- 
tained from 

N v = ($„|Sj[e Tt e T ] + Sl[e T "e T ]S v + [e Tt e T ]^|$„) 

= (* v \Sfa; + SfcS v +is;Sl\* v ) (3.17) 

with 



O = (e Ti Oe T ) f . c . + {e T 'Oe T ) . b . + (e T 'Oe T ) t . b . + 



and 



(3.18) 



7^= (e Tt e T )/. c . + (e I "e r ) . b . + (e'V'Jt.b. + ... (3.19) 

The f.c, o.b., t.b.,..etc abbreviations are used for the 
fully contracted, effective one-body, effective two-body 
...etc terms respectively Q. Terms containing only up to 
effective three-body diagrams will contribute to both the 
numerator and the denominator. The fully contracted 
terms are excluded on the basis of the linked-diagram 
theorem pj in the evaluation of the O and N. All the 
one-body terms have been taken into account as their 
contribution to the correlation effects is the largest. The 
dominant parts of the two-body terms have also been 
computed 0. Finally, these terms are contracted with 
Si and S v operators. 

Contributions from the normalization factor have been 
determined in the following way 



Norm = (V v \0\y v ){ 



1 



l + N v 



1}. (3.20) 



IV. RESULTS AND DISCUSSION 

We have used Gaussian type orbitals (GTOs) for the 
construction of single particle orbitals of the Dirac-Fock 
wavefunction(|$}), whose expression is given by |l9j | 



F; 



(L/S) 



(r) 



(4-1) 
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TABLE I: Magnetic dipole (A) and electric quadrupole (B) hyperfine structure constants in MHz of ground state and first 
excited state for Sc and Y 



States 


Scandium 
Theory (Ours) 

'a b 


Experiment \jx,6} 
A 


B 


States 


Yttrium 
Theory (Ours) 

A B 


Experiment [7. 8. 9. 10. Ill 

A~ B 


Sd"D 3/2 
3d 2 D 5/2 


268.67 27.45 
110.75 38.56 


269.556(1) 
109.032(1) 


26.346(4) 
37.31(10) 


4d 2 L> 5 /2 


57.34 23.37 
29.22 31.75 


57.217(15) 
28.749(30) 



TABLE II: Contributions from different coupled-cluster terms 
to the magnetic dipole hyperfine structure constant (A) in 
MHz, where cc stands for the complex conjugate part of the 



corresponding terms 




Scandium 


Yttrium 


CCSD(T) 
terms 


3d 2 D 3/2 3d 2 D 5/2 


4d 2 D 3/2 4d 2 D 5/2 


O (DF) 
O-O 
OSiv + cc 

OS2v + CC 

S\ v OSi v 
S\ v OS 2v + cc 

S\ v OS 2v + CC 


235.93 100.37 
4.82 2.41 

40.02 15.64 

-20.09 -22.93 
1.85 0.67 
1.21 -2.30 

25.79 25.59 


56.06 23.24 
2.37 1.04 
7.62 3.11 
-8.27 -5.08 
0.26 0.11 
0.32 -0.54 
3.62 9.45 


Imp 


nrtant effective two-body terms of O 


Sl v OTi + cc 
S\ v OT 2 + cc 
Norm. 


-9.07 -3.77 
0.00 0.00 
-10.94 -4.61 


-2.05 -0.85 
-0.32 -0.02 
-2.38 -1.19 



TABLE III: Contributions from different coupled-cluster 
terms to the magnetic dipole hyperfine structure constant (B) 
in MHz, where cc stands for the complex conjugate part of 
the corresponding terms 



Total 



268.67 



110.75 



57.34 



29.22 





Scandium 


Yttrium 


CCSD(T) 
terms 


3d 2 D 3/2 3d 2 D 5/2 


U 2 D 3/2 4d 2 D 5/2 


O (DF) 
O-O 
OSi v + cc 
OS 2v + cc 
S\ v OSi v 
Sl v OS 2v + cc 
S\ v OS 2v + cc 


23.34 32.76 
0.78 0.84 
4.05 5.53 
-0.97 -1.18 
0.18 0.24 
-0.24 -0.33 
2.74 3.99 


16.01 21.78 
0.69 0.76 
2.18 2.94 
3.90 5.78 
0.07 0.11 
0.06 0.99 
2.03 2.44 


Imp 


srtant effective two-body terms of O 


Sl v OTi + cc 
Sl v OT 2 + cc 
Norm. 


-0.91 -1.16 
0.00 0.00 
-1.14 -6.7 


-0.59 -0.79 
-0.01 -0.10 
-0.97 -1.30 


Total 


27.45 38.56 | 23.37 31.75 



with k— 0,1,2,... for s, p, d,- ■ ■ respectively. The function 
F^ S \r) stands for the large (L) and small (S) compo- 
nents of the Dirac wavefunction. c\ L is the expansion 
coefficient of the corresponding large and small compo- 
nents respectively. The kinetic balance condition pfjj has 
been imposed between the large and small components 
of the GTOs. For the exponents, the even tempering 
condition [21j 

Oi = 0^-1/3, i=l,-»,JV (4.2) 

has been applied. Here, N stands for the total number of 
basis functions for a specific symmetry. All orbitals are 
generated on a grid using a two-parameter Fermi nuclear 
distribution approximation given by 



where po is the average nuclear density, the parameter 
'c' is the half-charge radius, and 'a' is related to the skin- 
thickness which is defined as the interval of the nuclear 
thickness which the nuclear charge density falls from 



near one to near zero. 

In the present calculation, we have used a = 0.00525 
and (3 = 2.73 for all symmetries in both the systems. We 
have considered 35si/2, 30pi/2, 30^3/2, 30d 3 / 2 , 30d5/2, 
2O/5/2 and 20/ 7/2 and 38s 1/2 , 35p 1/2 , 35p 3/2 , 35d 3/2 , 
35^5/2, 2O/5/2 and 2O/7/2 GTOs to obtain Dirac- Fock 
wavefunction in Sc and Y respectively. All core electrons 
have been excited in the present calculations. 

We have calculated both A and B constants for 
3ci 2 £>3/ 2 and 3d 2 D 5 / 2 states for Sc and 4d 2 Z? 3 / 2 and 
4d 2 D 5 / 2 states for Y. We have used gi = 1.359 and 
Q = 0.22 for Sc and gi = 0.3268 and Q = 0.125 for Y 
respectively. All calculated results with corresponding 
experimental results are presented in table I. 

In table II and table III, we present contributions from 
different RCC terms for A and B values respectively. The 
first term O represents the DF value and is given in FIG. 
1 (i). The second term O — O is the contribution from 
the correlation between the core electros. OS\ v (FIG. 
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S 2 v f S 2v 

(ix) (x) 

FIG. 1: Goldstone diagrams representing hyperfine structure 
calculations 



tems with single s and p valence electrons. One of the 
strengths of coupled-cluster theory is to account for such 
intricate correlation effects to all orders in the residual 
Coulomb interaction. The two-body O and normaliza- 
tion terms make non- negligible contributions. 



V. CONCLUSION 

We have performed calculations of the hyperfine con- 
stants A and B for Sc and Y using RCC theory and 
high lighted the importance of electron correlation. Our 
results are in good agreement with experiments. This 
indeed demonstrates the power of RCC theory to cap- 
ture the interplay of relativistic and correlation effects 
for atomic systems with single d valence electrons. 



1 (ii)+(iii)) and OS 2v (FIG. 1 ((iv)+(v)+(vi)+(vii))) 
which represent pair-correlation and core-polarization ef- 
fects respectively make important contributions for both 
the atoms. It is interesting to note that the contributions 
of S f 2v OS 2v (FIG. 1 ((viii)+(ix)+(x)+(xi)+(xii) + (xiii))) 
is far more significant for Sc and Y than atomic sys- 
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